Dynamical Ordering and Directional Locking For Particles Moving Over 

Quasicrystalline Substrates 
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We use molecular dynamics simulations to study the driven phases of particles such as vortices 
or colloids moving over a decagonal quasiperiodic substrate. In the regime where the pinned states 
have quasicrystalline ordering, the driven phases can order into moving square or smectic states, 
or into states with aligned rows of both square and triangular tiling which we term dynamically 
induced Archimedean-like tiling. We show that when the angle of the drive is varied with respect 
to the substrate, directional locking effects occur where the particle motion locks to certain angles. 
It is at these locking angles that the dynamically induced Archimedean tiling appears. We also 
demonstrate that the different dynamical orderings and locking phases show pronounced changes as 
a function of filling fraction. 

PACS numbers: 05.60.Cd,74.25.Wx,05.45.-a,82.70.Dd 



Quasicrystals have received intense study since their 
discovery due to their unusual property of combining 
nonperiodicity with long-range order Recently there 
has been growing interest in understanding how interact- 
ingparticles such as vortices in type-II superconductors 
043] or charged colloidal particles [B|-[9| order in the pres- 
ence of a quasicrystalline substrate created with nano- 
lithographic or optical techniques. On a fully periodic 
substrate, these types of particles form states that are 
commensurate with the periodicity of the substrate at 
fillings which meet integer or fractional matching condi- 
tions flol - fl3j | . In the case of vortices, the commensurate 
fillings are associated with peaks in the critical current 
or the external force required to depin the vortices from 
the substrate 1(J 11 1. When vortices interact with qua- 



sicrystalline Penrose or decagonal pinning site arrays, a 
new set of peaks in the critical current appear at non- 
rational fields in addition to the commensurate fields due 
to a novel type of ordering of the vortices on these sub- 
strates @, IH . For colloids interacting with decagonal sub- 
strates, experiments show that strong substrates produce 
quasicrystalline order of the colloids but weak substrates 
permit the triangular lattice favored by the colloid-colloid 
interactions to form. Remarkably, for intermediate sub- 
strate strengths a new type of ordering arises consisting 
of a combination of triangular and square order. This 
has been termed an Archimedean- like tiling [5, |8fl . 

Numerous experimental and simulation studies have 
been performed on the dynamic phases of a collection 
of particles driven over a periodic substrate. One of 
the most striking observations is directional locking in 
which the particles cease to follow the direction of the 
drive and instead lock to a symmetry direction of the 
substrate. The locking effect produces a series of pro- 
nounced steps in the velocity versus drive angle curve 
l4 1 1 51 ] . Directional locking effects were experimentally 
observed for colloids moving through a periodic optical 
trap array 1(1 17 1 and for vortices driven at differing an- 



gles with respect to a pinning substrate [18] ■ For colloidal 
systems, directional locking can be used for practical ap- 
plications such as the fractionation of different species of 
particles, only one of which locks to the symmetry di- 



rection of the substrate [17J, [19|. An open question is 



whether directional locking can also occur for particles 
moving over quasiperiodic substrates. It has already been 
shown that the nonequilibrium, driven state for particles 
driven over random substrates can show dynamically in- 
duced ordered states which are inaccessible to the static 
system. For example, on strong random substrates which 
produce a disordered pinned state, the driven state can 
partially order to a moving smectic or to a triangular 
lattice @,[Il|. 



Here we show that directional locking and dynamical 
ordering can occur for particles moving over quasicrys- 
talline substrates. The dynamical ordering occurs only 
for certain driving directions, and changes the system 
from a quasicrystalline pinned state to a moving state 
with novel dynamical ordering that is similar to a square 
and triangular tiling. We call this ordering a dynamically 
induced Archimedean tiling. Other types of ordering to 
moving smectic, moving square, and moving liquid states 
are also possible as a function of filling and drive direc- 
tion. We specifically demonstrate the locking and order- 
ing for vortices in type-II superconductors, and note that 
the same results appear for charged colloids moving over 
decagonal substrates [2^ |. 

We model a two-dimensional system of N v interact- 
ing particles in the presence of a pinning array with 5- 
fold Penrose ordering, as shown in Fig. 1(a). In the 
superconducting system we consider a sample of size 
24A x 24A with periodic boundary conditions, where 
A is the London penetration depth. The vortex mo- 
tion evolves according to the following overdamped equa- 
tion: = F™ + Ff + Ff xt . Here R, is the lo- 
cation of vortex i and r\ is the damping coefficient. 
The repulsive vortex-vortex interaction force is F™ = 
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FIG. 1: (a) The locations of the Penrose tiled pinning sites 
(open circles) for a portion of the sample. Dashed lines in- 
dicate the positions of the Penrose tiles used to define the 
pinning locations, (b) The location of the vortices (filled cir- 
cles) and pinning sites (open circles) at BjB^, = 1.61 in the 
same sample at zero applied drive. Here £?0 is the field at 
which the number of vortices equals the number of pinning 
sites. Inset: The structure factor S(k) of the vortex positions 
from (b) has tenfold peaks indicative of quasicrystalline order. 




Ylf^i foKi(Rij / X)fiij ;, where K\ is the modified Bessel 
function, /o = 0q/(27t/ioA 3 ), 4>q — h/2e is the elementary 
flux quantum, Rij = |R,.j — Rj|, and R = (R^ — Rj)/ Rij. 
The substrate force from N p parabolic traps of radius 
r p = 0.35A and maximum strength F p ~ 1.85 has the 

form Ff" = fo(F p /r p )e(l - R lk /R p )R lkl where 6 

is the Heaviside step function. At the matching field 
of £?0, N v /N p — 1. Ff xt is the driving force from 
an external current applied uniformly to all particles. 
We initialize the particle positions using simulated an- 
nealing, then slowly turn on an external drive F ext = 
^(cos^x + sin^y) with 9 reported in degrees. We mea- 
sure the velocity in the y-direction V y = J2i=o v * ' V as 
we increment the drive angle 9. All forces and lengths 
are measured in units of /o and A. 

In Fig. 1(a) we plot the pinning site locations show- 
ing the 5-fold Penrose tiling and in Fig. 1(b) we show 
the positions of the non-driven pinned particle state at 
B/B<p — 1.61. The structure factor S(k) of the particle 
positions in this sample plotted in the inset of Fig. 1 (d) 
has a 10-fold peak structure, indicative of quasicrystalline 
order such as that found in the strong substrate limit for 
colloidal systems In general, for this F p the pinned 
states at different filling fractions have quasicrystalline 
ordering. 

Figure 2(a) shows (V y ) versus the drive angle 9 for a 
sample with B/B^ = 3.225 and F D = 2.0. At this drive, 
all the particles are in motion. The clear set of steps in 
(V^) is a signature of the directional locking that occurs 
when the particles move alon g a single direction over a 
range of drive angles [I3 - fl7l . Ixol] . For the same set of 
parameters but with randomly placed pinning, no steps 
occur in (V y ) (22}. The locking angles are related to the 
tenfold orientational ordering, and occur at multiples of 



FIG. 2: (a) The average velocity in the y-direction (Vy) vs 
drive angle 6 for the system in Fig. 1(a) with B/B$ = 3.225 
at Fd = 2.0. A series of pronounced steps appear at multiples 
of 36°, marked 0/1, 1/1, and 2/1, when the particle motion 
becomes locked to certain symmetry directions of the sub- 
strate. Locking at fractional multiples of 36° produces many 
smaller steps. Inset: A blow up of the main panel near the 
1/1 step showing some of the fractional steps, (b) The corre- 
sponding Pa vs 6 showing that the system is more ordered on 
the locking steps. 

9 = 360°/10 = 36°, as highlighted in Fig. 2(a) for 9 = 0° 
(the 0/1 step), 36° (the 1/1 step), and 72° (the 2/1 step). 
There are numerous smaller locking steps associated with 
rational fractions of 36°, including 1/4, 1/2, 2/3, 3/4, 
5/4, and 5/2. Some of these smaller steps are visible in 
the inset of Fig. 2(a), which shows a blowup of the main 
panel. This result proves that it is possible for directional 
locking to occur on quasiperiodic substrates, indicating 
that orientational order rather than translational order 
of the substrate is the essential ingredient for directional 
locking. This opens the possibility of creating novel frac- 
tionation devices using quasiperiodic substrates, and may 
also be relevant to understanding frictional studies per- 
formed with quasicrystalline substrates. We find that on 
the locking steps, there is a higher amount of order in 
the particle lattice. In Fig. 2(b) we plot the fraction of 
sixfold coordinated particles Pe = J2i " $( z i — 6) as a 
function of 9, where the particle coordination number Zi 
is obtained from a Voronoi construction that is not sensi- 
tive to square orderings. In the pinned state Pq = 0.525, 
while for the driven systems Pe reaches its highest val- 
ues on the locking steps and is lower in the non-locking 
ranges of 9. 

We tessellate the positions of the particles in the mov- 
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FIG. 3: Delaunay triangulations of the particle positions 
in the sliding state. Bonds longer than l.lao are omitted, 
where ao is the lattice spacing of a triangular lattice with the 
same density. Yellow (light) polygons are triangular and blue 
(dark) polygons are square or nontriangular. For a sample 
with B/B^ = 3.225 and F D = 2.0 we show (a) the 0/1 step 
at 9 = exhibiting the oriented triangles and squares of the 
dynamically induced Archimedean tiling, (b) the 1/1 step at 
9 = 32° where the Archimedean ordering aligns with the di- 
rection of motion, and (c) 9 = 41° in the unlocked region just 
above the 1/1 step, (d) The triangulation for the 1/1 step at 
9 = 32° for a sample with B/B^, = 1.61 and Fd = 2.0, where 
most of the tiles are square. Smectic-like ordering appears in 
S(k) corresponding to the Archimedean tilings of (e) panel 
(a) and (f) panel (b). (g) S(k) for the unlocked state shown 
in panel (c) has a ring structure indicating liquid ordering, 
(h) At the lower field of B/Bcj, = 1.61, S(k) corresponding to 
panel (d) has square ordering. 



ing state as in Refs. 0, Q in order to reveal the nature of 
the dynamical ordering. In Fig. 3(a) we show the tessel- 
lation of the moving state on the 0/1 step from Fig. 2(a), 
where the flow is locked in the ^-direction. Both square 
and triangular tiles are present, and the square tiles are 
aligned with the flow along the x direction. This state is 
very similar to the Archimedean type tiling found for col- 
loids pinned by substrates of intermediate strength 0, [|| . 
The smectic features of the corresponding structure fac- 
tor S(k) in Fig. 3(e) indicate that unlike the static col- 
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FIG. 4: (a) The depinning threshold F c Tr indicating the 
width of the 0/1 locking step vs B/B^. (b) Pa vs B/B^. 
The broad peak in both quantities for 1.7 < B/B,p < 4.3 
appears when the system forms the dynamical Archimedean 
tiling illustrated in Fig. 3(a), while the dip at B/B^ = 1.6 cor- 
responds to the moving square structure shown in Fig. 3(h). 
Additional peaks in P§ occur at B/B^ = 1.0 and B/B^ = 
1/r, where r is the golden mean. Dashed lines and the labels 
a-d indicate the fields shown in Fig. 3(a-d). 



loidal tiling, the dynamical tiling does not have a one- 
dimensional (ID) quasiperiodic structure, which would 
appear as a series of peaks in S(k) [||. More recent colloid 
experiments reveal the development of smeared smectic- 
like peaks in S(k) similar to those in Fig. 3(e) when the 
strength of the laser-induced substrate is increased 
Moving smectic states also appear in systems with ran- 
dom pinning, but these states have a purely triangular 
tiling containing some dislocations dHHH]. 

We find that oriented Archimedean tilings appear on 
each locking step, as illustrated for the 1/1 step in 
Fig. 3(b) and (f). T he p articles move along ID channels 
on the locked steps [22|, but follow winding and rapidly 
changing trajectories when not on the steps. In Fig. 3(c) 
we show the tessellation for driving just above the 1/1 
step in a disordered flow regime. The orientational or- 
dering is lost and the corresponding S(k) in Fig. 3(g ) has 
a ring structure characteristic of a moving liquid [20j . We 
find that the relative number of square and rectangular 
tiles appearing on the locking steps varies as a function 
of filling. For example, at B/B^ = 1.61 on the 1/1 lock- 
ing step, Fig. 3(d) shows that the tiling is predominantly 
square, and S(k) in Fig. 3(h) has square ordering. 

To understand where the different moving phases oc- 
cur as well as the changes in the effectiveness of the di- 
rectional locking, in Fig. 4(a) we plot the transverse de- 
pinning force F^ r determined by F^ r = Fd sin# e , where 
9 e is the end of the 0/1 locking step, versus B/B^ and 
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in Fig. 4(b) we show the corresponding Pq versus B/B^. 
There is a broad maximum in F^ r associated with the 
Archimedean type ordering for 1.7 < B/B^ < 4.3 and 
another sharper peak in F*£ T at B/B<p = 0.62. Both 
of these features are accompanied by peaks in Pq. Our 
transverse depinning measurements, in which the vortices 
are depinning from a moving state, differ in several ways 
from the longitudinal depinning from a static state previ- 
ously measured in simulations and experiments @, 0, [23| ■ 
For example, Fj r shows no peaks at B/B^ = 0.81 or 1.0, 
although we do hnd a weak peak in in Pg at B/B^ = 1.0. 
In agreement with the experiment of Ref. Q , we observe 
a prominent peak in F^F r at B/B^ = 0.62 = 1/r where 
t = (1 + v / 5)/2 is the golden mean. The lowest value 
of Fj r occurs at B/B^ = 1.61 w r. At this field, the 
widths of the locking steps are suppressed and the sys- 
tem has the square ordering shown in Fig. 3(d,h). For 
B/Brj, > 4.0, the ordering of the locked phases is reduced 
and S(k) develops a ring structure on the locking steps, 
while the vortex- vortex interactions are strong enough in 
the nonlocking regimes to produce sixfold ordering. For 
0.5 < B/B<f, < 1.5 the system generally exhibits a smec- 
tic ordering containing few square tiles corresponding to 
a highly defected triangular lattice similar to that found 
for driving over random substrates. 

We find that the ordered moving phases only occur 
when Fd is sufficiently large. For Fd < 1-25 with 
F p = 1.85, significant plastic flow occurs, S(k) is al- 
ways liquidlike, and the step structures in the velocity 
curves disappear. For Fd F p the width of the steps 
slowly decreases with increasing Fd- We also find the 
same results reported here for driven colloids interacting 
with a screened Coulomb interaction and a quasiperiodic 
substrate 22j . Sevenfold or tetradecagonal quasiperiodic 
substrates also produce directional locking, but the lock- 
ing steps are significantly reduced in width [22 1. 

We have shown that interacting particles such as vor- 
tices or colloids driven over a quasiperiodic substrate ex- 
hibit rich dynamical behaviors including pronounced di- 
rectional locking where the particles prefer to move at 
particular angles with respect to the substrate. Direc- 
tional locking has already been studied in periodic sys- 
tems, but this study is the first demonstration of direc- 
tional locking on quasiperiodic substrates. For strong 
substrates, the pinned state has quasicrystalline ordering, 
but the moving state can have square or smectic order- 
ing depending on the orientation of the drive. At certain 
filling fractions the system forms a moving Archimedean 
tiling ordered state similar to the pinned Archimedean 
tiling ordering observed for colloids on decagonal sub- 
strates, but with smectic rather than quasiperiodic char- 
acter. We also show that the dynamically ordered states 
produce distinct signatures in the transverse depinning 
threshold. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 



DE-AC52-06NA25396. 
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SUPPLEMENTARY MATERIAL 



Directional Locking for Vortices on Tetradecagonal 
Quasicrystalline Substrates 

Here, we show that the directional locking effects de- 
scribed in the main text for driven vortices moving over 
decagonal quasicrystalline substrates also occur when 
the vortices are driven over tetradecagonal or sevenfold- 
symmetric quasicrystalline substrates such as that shown 
in Suppl. Fig. 1. We use the same pinning parame- 
ters as for the decagonal substrate with F p — 1.85 and 
r p = 0.35A, and apply a drive Fd = 2.0 at B/B^ = 3.9. 
In Suppl. Fig. 2(a) we plot (V y ) versus the drive angle 
9 and in Suppl. Fig. 2(b) we show the fraction of six- 
fold coordinated particles Pq versus 9. There is a clear 
set of steps in (V^) which are accompanied by enhanced 
sixfold ordering as shown by the increases in Pq. For 
the tetradecagonal substrate potential, the major lock- 
ing steps occur at drive angles that are integer multiples 
of 360°/14, producing the 1/1 step at 25.7°, the 2/1 step 
at 51.42°, and the 3/1 step at 77°. For the decagonal 
substrate the locking steps fall at integer multiples of 
360°/10, as described in the main text. Fractional lock- 
ing steps are nearly absent for the tetradecagonal sub- 
strate; however, for lower fillings it is possible to resolve 
some fractional steps. This is shown in Suppl. Fig. 3(a,b) 
for B/B^ = 2.96. The fractional steps m/4 with m in- 
teger are the most pronounced and are associated with 
enhanced sixfold ordering in Pq. This result indicates 
that directional locking is a generic feature for parti- 
cles moving over different types of quasicrystalline sub- 
strates. The width of the locking steps is generally much 
narrower for the tetradecagonal substrate than for the 
decagonal substrate. Additionally, although the decago- 
nal substrate produced pronounced five-fold ordering of 
the particles for some fillings, we do not find strong seven- 
fold particle ordering on the tetradecagonal substrate. 
Instead, on the locking steps, we observe smectic ordering 
composed of sixfold-coordinated particles and a limited 
number of dislocations. 



Directional Locking for Colloids on Decagonal and 
Tetradecagonal Substrates 
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FIG. 1: The locations of the pinning sites (open circles) 
for a portion of the sample with seven-fold or tetradecagonal 
ordering. Dashed lines indicate the positions of the tiles used 
to define the pinning locations. 




We next show that the same directional locking effects 
described in the main text for vortices driven over qua- 
sicrystalline substrates also occurs for colloidal particles 
driven over quasicrystalline substrates. We model a two- 
dimensional system of N c interacting colloids in the pres- 
ence of fivefold and sevenfold quasicrystalline substrates. 
The substrate is composed of localized pinning traps with 
maximum strength F p and radius r p = 0.35. We simulate 



FIG. 2: (a) The average velocity in the y-direction (Vy) vs 
drive angle 8 for the system in Suppl. Fig. 1 with a seven- 
fold quasicrystalline pinning array for F p = 1.85, r p — 0.35A, 
B/ =3.9, and Fd = 0.2. Several steps appear at the direc- 
tional locking angles which are integer multiples of 360°/ 14. 
The 1/1, 2/1, 3/1, and 4/1 lockings are clearly visible, (b) 
The corresponding Pa vs 6 shows that on the locking steps the 
system develops a considerable amount of sixfold ordering. 




e 

FIG. 3: (a) (V y ) vs drive angle 6 for the system in Suppl. 
Fig. 2 but for a filling of B/Btt, = 2.96 where additional frac- 
tional locking steps are visible, particularly at fractions of 
m/4 where m is an integer. Some representative steps are 
marked, (b) The corresponding Pa vs 0. 



the motion of the colloids using the same procedure used 
previously to model colloid dynamics on random sub- 
strates [ij and periodic substrates [3] , by integrating the 
following equation of motion: rjdRd/dt = F^ c + Ff + Fd- 
Here Ri is the location of colloid i and rj is the damp- 
ing coefficient. The colloid-colloid interaction potential 
has the Yukawa form V(Rij) — (E /Rij)exp(—KRij), 
where Rij = |Rj - TLj\, E = Z* 2 /(4iree ), e is the 
solvent dielectric constant, Z* is the effective charge 
of each colloid, and 1/k is the screening length. The 
colloid- colloid interactions are repulsive and are given by 
Ff c = — Y^fli VV(iJy). Lengths are measured in units 
of ciq and forces in units of Fq = Eo/ao. The substrate 
force term Ff comes from N p pinning sites placed in a 
decagonal or tetradecagonal pattern and has the same 
form as described for the vortex system. The colloid den- 
sity relative to the pinning density is N c /N p . We neglect 
hydrodynamic interactions in the colloidal system and 
assume that the colloids are in the low volume fraction, 
highly charged, electrophoretically driven limit. The ex- 
ternal force Fd is the same as that used for the vortex 
case and we measure (V y ) and Pq as a function of the 
drive angle 9. 

In Suppl. Fig. 4(a) we plot (V y ) versus 9 for a col- 
loidal system on a decagonal substrate with F p = 0.75 
at N c /N p = 2.9, and in Fig. 4(b) we show the corre- 
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FIG. 4: The directional locking for charged colloidal particles 
driven over a decagonal substrate in the same manner as the 
vortices in the main text. Here N c /N p = 2.9 and F p = 0.75. 
In general, the same features observed for the vortex sys- 
tem also appear for the colloidal system. (a){V y ) vs 9 with 
the main locking steps at integer multiples of 360°/ 10 high- 
lighted at 1/1 and 2/1. Several fractional steps also appear 
at rational fractions of the integer steps such as at 1/4, 3/4, 
and 3/2. (b) The corresponding P§ vs 9 shows that along the 
steps the system has higher ordering. On the main integer 
matching steps, the system forms the dynamically induced 
Archimedean ordering also observed in the vortex system (as 
described in the main text). 



sponding Pq versus 9. Here, all the general features of 
the directional locking in the vortex system also appear 
for the colloidal system, including the dominant locking 
steps at integer multiples of 360°/10 such as 1/1 and 2/1 
and additional fractional lockings at 1/2, 3/2, 5/2, 1/4, 
and 3/4. Each locking step is associated with enhanced 
ordering of the colloidal lattice when the system forms a 
dynamically ordered Archimedean tiling state. 

In Suppl. Fig. 5(a,b) we show (V y ) and P§ versus 9 for 
colloids moving over a tetradecagonal substrate. We find 
the same features illustrated in Suppl. Fig. 3 for vortices 
moving over a tetradecagonal substrate. The dominant 
locking steps occur at integer multiples of 360°/14. In 
general, the locking steps are weaker than for the decago- 
nal substrate and the dominant fractional steps occur at 
ratios of m/4 with m integer. The fractional steps are 
also associated with a series of peaks in Pg, where addi- 
tional fractions beyond m/4 appear as smaller peaks that 
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do not show full ordering. These results indicate that 
the directional locking effects on quasicrystalline periodic 
substrates can be observed for various types of interact- 
ing particles including colloids and vortices, and that the 
directional locking is a robust feature in these systems. 
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FIG. 5: Directional locking for charged colloidal particles 
driven over a sevenfold or tetradecagonal substrate for the 
same parameters in Suppl. Fig. 4. (a) (V y ) vs 6 and (b) 
P(, vs 8. The same features found for the vortex system in 
Suppl. Fig. 3 appear for the colloid system, including locking 
steps falling at integer multiples of 360°/ 14. We also observe 
the strongest fractional locking steps at fractions m/4 with m 
integer. 



